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ABSTRACT 


This paper extends an analysis of plastic stress waves, originated hy G. 
I. Taylor in reference 1, for Gylindrical metallic projectile in impact to an 
analysis of a hemispherical shell suffering plastic deformation during the 
process of impact. In that, it is assumed that the heMspherical shell with a 
prescribed launch velocity impinges a fixed rigid sphere of diameter equal to 
the internal diameter of the shell. Particiilarly this- study is directed in 
order to investigate the dynamic biaxial state of stress present in the shell 
during deformation. 


The results of this analysis are compared with Taylor's reference 1 and it 
has been found that this analysis is an extension of the one-dimensional 
analyses of references 1, 2, 3, and 4, to spherical coordinates^ It is valuable 
for studying the state of stress .during large plastic deformation of a hemi- 
spherical shell. 


INTRODUCTION 


The object of this paper is to develop an analysis of plastic hemispheri- 
cal stress-wave propagation and to use this analysis for determining the 
dynamic biaxial yield stress, The.Tresca yield, criteria is used as the yield 
condition. Higher order terms are included in the derivations; thus, this 
analysis is valid for large deformations. 

G. I. Taylor in reference 1 used the governing physical laws and the 
geometry during plastic Reformation of the cylindrical projectile to formulate 
differential equations which are solved in order to determine the dynamic yield 
stress in impact. This analysis of a hemispherical shell impacting a fixed 
rigid sphere, of diameter equal to the internal diameter of the shell, is 
similar to the analysis of a cylindrical projectile impacting a rigid target of 
references 1, 2, 3, and 4. In fact, in all these cases during impact it is 
assumed that when the stress rise exceeds the elastic limit of the material. 
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two waves are generated. The first of these is the elastic wave, which travels 
with a velocity c. It is followed by the plastic stress wave which travels 
with a slower velocity v. Through an analysis of the propagation of these two 
stress waves, a method is formed which can be used to determine the dynamic 
yield stress of the material of a hemispherical shell. The proper choice of 
the time increment, dt, simplifies the analysis greatly. The choice is to make 
the time increment eq.ual to the length of time required for the elastic wave to 
complete a double passage of the elastic zone. If the difference equations are 
derived by utilizing this time increment, which is eliminated by combining the 
derived difference equations, the governing equations which are derived are 
free of this time increment. This mathematical approach, for the biaxial state 
of stress of the hemispherical shell, closely parallels Taylor's analysis of 
the cylindrical projectile. 


nomenclature 

Projection at the elastic-plastic boundary undeformed area at time t 
A Projection at the elastic-plastic boimdary deformed area at time t + dt 
a Initial inside radius 

b Initial outside radius 

c Elastic wave Velocity 

dh Incremental plastic radius 

dr Incremental elastic radius 

dt Time for a double passage of elastic region by elastic wave 
E Young's modulus 

h Thickness of plastic region at time t 

r^ Initial elastic length of shell in the radial direction = b - a 
r^ Final total length of shell in the radial direction 

r Thickness of shell in the elastic zone at time t 

R Final thickness of the elastic region 

Dyhamlc yield stress - calculated by the approximate method 
S Dynamic yield stress - calculated by the exact method 
t Time 
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U Initial radial velocity due to launch velocity 
u Particle velocity in the elastic region 

V Absolute velocity of the plastic wave front 

V Yield stress in uniaxial tension 

e^ Initial radial strain 

e Radial strain at time t 

p Density 

V Poisson’s ratio 


ANALYSIS 

Problem Description 

A hemispherical metallic shell strikes with a prescribed velocity a rigid 
sphere (of diameter equal to the internal diameter of the shell) which is 
permanently fixed at a base retaining zero velocity during the process of 
impact. During this impact, a radial motion is directed from the internal sur- 
face of the shell toward the external surface of the same( The radial particle 
velocity of the internal surface of the shell is initially the same as the 
impact velocity and is denoted by U. If the biaxial stress exceeds the elastic 
limit, two waves are generated at the internal surface of the shell. The first 
wave is the elastic wave which travels with velocity c. The second is the 
plastic wave which travels with velocity v.. The elastic compressive stress 
wave, which propagates radially outward in the elastic region with velocity c, 
will reduce the impact velocity U to U-(S/pc). During this time, the stress 
reaches the elastic limit. This elastic wave will reflect at the external 
surface of the sphere, resulting in an elastic tensile wave being superposed on 
the compressive elastic wave. The material which has been passed by this re- 
flected elastic wave is stress free and has a velocity equal to U-(2S/pc). At 
the particular time when this wave reaches the elastic-plastic boundary, the 
shell is in a condition similar to the initial impact, except that its speed 
is equal to U-(2S/pc) and its elastic thickness is less than the original value. 
At this time, it is assumed that the plastically deformed material will be 
attached to the sphere and acts on the elastic part of the shell as a rigid 
material. This continues until the speed of the shell becomes equal to zero. 
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Assumptions 


In order to work out the mathematical analysis of this prohlem, several 
"basic assumptions are needed. 

First, for axially symmetrical analysis, the shell must "be ssnnmetric with 
respect to its axis of symmetry and maintain this symmetry during the process 
of impact. The second assumption is that the elastic strain is negligi"ble. 

This asstimption is valid if the plastic strain is large, thus nmiking the elas- 
tic strain very small in comparison with the plastic strain. Along the same 
lines as the previous assumption, the third assumption is that the material is 
taken to "be perfectly plastic. Although no material "behaves exactly in a 
perfectly plastic manner, some materials approach this type of behavior at high 
strain rates. This dynamic-plastic stress-wave analysis is for extremely high 
strain rates. Thus it is possible to assume that the material is perfectly 
plastic without the loss of much generality in the solution. The fourth 
assumption, which is usually made in plasticity problems, is that the density 
of the shell material remains constant. The fifth and final assumption is 
that the material in the plastic region, after being deformed, does not possess 
elastic properties; thus it behaves as a rigid material with zero velocity. 


Physical laws 


By considering the problem description, and assumptions, the governing 
physical laws can be formulated. 


Choose the time increment, dt, to be equal to the time required for a 
complete double passage of the elastic wave through the elastic region. Since 
the length of material in the elastic zone is defined as r and the elastic wave 
velocity is c, it follows that 

dt = 2r/c (1) 


where 


c = {E (1 - v)/ {p(l + v)(l - 2v)}} 
dh = V (2r/c) 
dr = -(u + v) (2r/c) 
du = - 23/ (pc) 


1/2 


( 2 ) 

(3) 

(V) 


Using equation (l) to eliminate c from equations (2), (3)» and (4) results in 
dh 


dt 


= T 


(5) 


dt 


-(u + v) 


( 6 ) 
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f =-S/(pr) 

for conservation of mass 

A V = (u + v) A 
o 

The momentum equation reduces to 


(7) 

( 8 ) 


S (a - A ) = 1/2(A + A ) (u + v) u p 
o o 


(9) 


The radial strain is defined, at the plastic houndary, hy 

e = 1 - A /A 
o 

Combining equations (8), (9), and (lO) 
pu^/S = 2 e^(2 - e) 

Combining equations (6), (7), (8), and (lO) 
dr 


du 


= pur/(Se) 


Integrating equation (12) 
Log^(r^) 


■/ 


1/e d{e^/(l - e/2)} 


= i^/(2 - e) - 2 Log (l - e/2) + Constant 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


At time t = 0, u = U, r = r , and e = e j thus equation (ll) and equation 
(13) become, respectively 


pU^/S = 2 e^ /(2 - e ) 
1 1 


(IV) 


Log^(r/r^)^ = k/{2 - e) - 2 Log^ (l - e/2) - k/{2 - e^) 
+ 2 Log^ (1 - e /2) 

e 1 

When all motion has ceased, r = R, and e = 0, and R can be measured. 

Log (R/r )2 = 2 - 4/(2 - e ) + 2 Log (l - e /2) 
e o 1 1 

Combining equations (5), (6), (8) and (lO) 

r 

h = I dh = - I (1 - e), dr 


/*-i 


(15) 

(16) 

(17) 
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Combining equations (7), (ll)j and (lU) 

Ut/r = e (l - e f ^r/r (l - gA)/(1 - z/2)^^ (l8) 

o 1 1 J ° 

£ 

If uniformly spaced values of are placed in equations (ill-) and (l6), 
pU^/S vs R/r can be plotted. Evaluation of equation (iT) for h is accomp- 
lished by Simpson’s rule integration. Results of these calculations are 
plotted in Figiire 1. 

Two different methods of integration were employed to evaluate the inte- 
gral equations (17) and (l8). The first method used was a Simpson's rule 
integration. Results for various are plotted in Figure 2, 

The second method of integration was using the asymptotic expansion of 
the integrals. References 5» 6, eind T provide information on asymptotic 
power-series expansions,. Values obtained by asymptotic expansion agreed well 
with those obtained by Simpson’s rule integration. 


To develop a simple formula for calculating the dynamic yield stress from 
measurements made before and after the impact, it will be additionally assumed 
that the plastic boundary propagates at a constant velocity from the inside 
radius a to its final position. The velocity of the plastic boundary equals C. 

Combining equations (6) and (j) 

= S/{pr(u + C)} (19) 

Integrating equation (19) results in 

S/p Log (r/r ) = l/2u^ + C u - 1/2U^ - C U (20) 

0 o 

When u = 0, r = R and equation (20) becomes 

S/p Log (R/r ) = - 1/2U^ - C U (21) 

0 0 

At time t = 0, u - U. Assuming u decreases to zero uniformly with time, 
in a time equal to T 

T = (r - R)/C = 2(r - r )/U 

1 0 1 

Rearranging C/U = 1/2 (r - R)/(r - r ) 

1 o 1 

Therefore, equation (21) becomes 

S^/pU2 = (r^ - R)/[2(r^ - r^) Log^ (r^/R)] (22) 

The fact that the decrease in u is not uniform results in an ehror which 
can be calculated. Combining equations (3) and (20) 
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= 2S/p Logg 

When all motion has ceased, u = 0, r = R, and e = 
■becomes 

2S/p Log (R/r ) = C2 - (U + C)^ 
e o 

Letting 

2S/p = a2 

K = (U + C)/a 

Ri = r/r 
^ o 

ti = at/r 

o 

Tj = aT/r^ 

where T is the time from the initial impact until 
equals zero 

= (k 2 + Log Ri)' 

ati ® 



(23) 

0. Therefore, equation (2l) 

(2U) 


the plastic zone velocity 


so that 


i 

1 = / „ 2 . *'^^6 "I’dR 

•tvT^r i^) - 


exp[(-) - K^3 

Gt 


- 1/2 
1 


Letting 


+ Logg Rl 


(25) 


results in 



(26) 


6 . 


Values of F(k) = e“' 
Equation (26) can be 


I 


2j2 

e dZ have been tabulated in references 5 and 
expanded using this function, P(K) , to give 


Ti = 2 {F(K) - Exp (-K2 + (C/a)2) F(C/a)} (27) 

Previously it was assumed that the plastic boundary moves with a constant 
velocity C, or 


C T = r^ - R 
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or in dimensionless form 


C/a Ti = ri/r - R/r (28) 

■ o o 

Rearranging, equation (2^1) becomes 

Logg (r^/R) = - (C/a)2 (29) 

Since R/r^ and r^/r can be measured. C/a, U/a, K and T can be evaluated 
from equations (27), (28)° (29), and 

K = U/a + G/a (30) 

Combining equations (27) s (28), and (29) 

r/r =2 C/a F(K) - {2 C/a F(C/a) - 1} R/r (3l) 

I o o 

Since 

2S/pU2 s cx2/u2 = 1/(K - C/a) 2 

dividing this equation by equation (22) therefore results in 

S/S = (r - r )/(r - R) [Log (r /R)/(K - C/a)2] (32) 

1 0 1 o e o 


Due to the complexity of these equations, the correction factor, S/Si, 
cannot be determined directly. To determine S/Si given R/r and r^/r^ it is 
easiest to first form the curves of S/Si vs r^/r 'with contours of equal h/r . 
Values for this curve can be obtained by taking a value of R/r and values o? 
C/a which cover the desired range. Therefore, using equation ?29), equations 
(3l) and (32) can be eval\xated. 

The asymptotic expansion of F(K) is 

F(K) = 1/(2K) + 1/(Uk3) + 3/(8k5) + 15/(i6k7) + .... (33) 

Through some complex manipulations, it can be shown, although it will not 
be presented here, that as C/a ->■ » 

r 1 /r -»■ 1 . 0 ( 3L ) 

i. o 

and 

S/Si = 2 {1/(1 - R/r^) - l/(Logg(r^/R))} (35) 

. - f 

From this equation, the limiting values of S/Sj can be determined as 

ri/r ->■ 1.0. 

^ o 

This completes the analysis of the problem. Thus, if values of r , r^, 
and h are given, the dynamic yield stress can be calculated for the hemispheri- 
cal shell . 
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DISCUSSION AND CONCLUDING EEMAEKB 


In this paper, a method is developed to investigate the propagation of 
plastic stress waves in a hemispherieal shell. In particular, this stu(3y 
investigates the dynamic yield stress due to the impulsive loading Initiated 
at the Interface of the shell. This mathematical approach, for determining 
the biaxial state of stress of the hemispherical shell, closely parallels 
Taylor's analysis of the cylindrical projectile. It is interesting to note 
that if all higher order terms were dropped from this analysis, the results 
would be the same as those of reference 1 except that the components are 
defined differently. Graphs which are drawn from this analysis in Figures 1 
and 2 are similar to Figures 2 and 3 of reference 1. In addition, comparison 
between the results of this analysis and the analysis of reference 1 is 
possible. In fact, when these two analyses are compared, one can observe that 
the results of the present work parallel the experimental data more closely, 
than the restilts of reference 1. This is due to the fact that one-dimensional 
analysis may not possibly explain the spreading out of the projectile near the 
target. This phenomenon requires taking into account the inertia in the radial 
direction. 


The derivation of the yield stress correction factor is almost identical 
with the results of reference 1 on page 29T* Singularities were observed 
which were not discussed in reference 1. The discontinuities occurred just 
before f^/r -> 1,0. If the discontinuity is Ignored, the results are similar 
to those or reference 1. 

A method has been presented by which the dynamic yield stress can be 
calculated, using the Tresca yield criteria, from the radial, expansion of a 
hemispherical shell. The approximate yield stress can be calculated from 
equation (22), if the initial conditions, final conditions, U, and p are 
specified. The dynamic yield stress could also be calculated from Figirre 1, 
Thus, the dynamic yield stress can be determined if certain initial and final 
experimental conditions are specified, including the launch velocity, density, 
and geometrical considerations of the shell. The motion of the plastiG 
boundary, as shown in Figure 2, is similar to the results obtained in Figure 4 
of reference 8. Their choice of coordinates is dlfferei^t, which accoimts for 
many of the differences between the shape of their curve and of Figure 2. 
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